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GENERAL  NOMENCLATURE  (Continued) 

2 2 

v , R ■ coefficient  of  determination 
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SECTION  1 

(This  Section  is  Unclassified) 

INTRODUCTION 

The  objective  was  to  construct  improved  height  of  burst  contours  for 
dynamic  pressure  impulse.  In  particular,  we  wished  to  extend  the  set  of  con- 
tours previously  generated  (Reference  1)  to  higher  values  of  dynamic  pressure 
impulse  and  to  obtain  data  points  along  the  contours  at  additional  scaled 
burst  heights,  i.e.,  to  better  determine  the  contour  shapes.  To  do  this  we 
analyzed  all  of  the  data  available.  This  includes  dynamic  pressure-rime  wave- 
forms for  shots  not  used  in  Reference  1 (because  no  trucks  were  exposed  on  these 
shots):  TEAPOT  Homet-S,  Post-11,  and  Zucchini-14;  PUJMBBOB  Franklin-2,  Wilson-4 
Hood-6,  Kepler-9,  and  0wens-10.  It  also  includes  data  for  shots  in  which  tanks 
and  howitzers  were  exposed  but  no  dynamic  pressure  measurements  were  made: 
GREENHOUSE  Easy- 2;  TUMBLER-SNAPPER  Fox-6  and  How-7;  UPSHOT -KNOTHOLE  Annie-1*, 
Nancy-2,  Badger-S,  and  Siaon-7. 


* There  was  one  dynamic  pressure  measurement  on  this  shot. 
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SECTION  2 

(This  Section  is  Unclassified) 

ANALYSIS  OF  BLAST  MEASUREMENTS 
2.1  Measurements  For  Evf.tf  Not  Previously  Considered. 

We  consider  first  the  shots  not  previously  considered  for  which  some 
dynamic  pressure  measurements  were  available.  Some  of  the  waveforms  were  not 
very  good  owing  to  noise  (mixed  inextricably  with  real  physical  irregularities), 
baseline  movement,  or  gage  record  cutoff.  However,  many  of  the  waveforms  were 
found  useable,  being  a quality  similar  to,  or  not  significantly  inferior  to, 
those  used  in  Reference  1.  General  information  on  all  of  the  nuclear  tests  of 
concern  in  this  report  is  given  in  Table  1.  Dynamic  pressure  impulse,  scaled 
and  unsealed,  is  given  in  Table  2 for  each  Operation/Event,  scaled  height  of 
burst,  and  ground  range  for  which  we  have  obtained  results.  Also  given  in 
Table  2 is  the  manner  in  which  we  obtained  the  values  of  dynamic  pressure 
Impulse . 

We  discuss  the  data  for  the  shots  in  the  order  listed  in  Section  1. 

TEAPOT,  Hornet-5.  As  indicated  in  Table  1 we  used  BRL  gage  data.* 
There  were  two  gage  results  at  each  of  the  three  ground  ranges.  In  two  instan- 
ces the  two  results  are  in  good  agreement.  At  the  closest-in  range,  256  meters, 
the  two  gages  differ  by  about  a factor  of  two  but  one  of  the  two  appears  to  be 
correct,  the  other  badly  in  error.  We  used  the  average  of  the  two  results  for 
the  329  and  460  metre  ranges,  and  the  apparently  correct  results  for  the  256 
metre  ground  range. 

TEAPOT,  Post-11.  For  this  shot  dynamic  pressure  waveforms  are  avail- 
able in  Reference  2.  We  used  a planimeter  to  obtain  the  area  under  each  curve 
and  then  used  the  data  reduction  procedure  described  in  detail  in  Reference  1. 

A detailed  discussion  of  accuracy  of  this  procedure  is  also  given  in  Reference 
1.  For  this  shot  waveforms  are  given  at  four  ground  ranges  in  Reference  2. 

The  waveforms  are  not  good  but  are  of  sufficient  quality  to  be  useful.  As 
indicated  in  Table  2 we  used  the  average  of  BRL's  result  and  our  result  for 
the  two  ground  ranges  at  which  BRL  provided  results.  Our  results  and  BRL's 
results  agreed  to  within  about  3%  which  is  quite  good;  (the  two  results  are 
based  on  the  same  gage  data  but  the  data  reduction  procedures  are  independent; 
we  have  analyzed  errors  inherent  in  data  reduction  in  Reference  1,  and,  in 
many  instances,  errors  can  be  quite  large). 

For  the  two  ground  ranges  at  which  BRL  did  not 
reduce  the  data,  we  used  the  result  of  our  own  data  reduction  process. 

TEAPOT,  Zucchini-14.  For  this  shot  we  averaged  our  result  with 
the  result  of  BRL's  data  reduction  process  for  the  three  listed  ground  ranges. 


* BRL  gage  data  used  in  this  report  are  published  in  DNA-TR-85-161 . 
Furnished  by  Mr.  J.  Keefer,  BRL. 
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In  each  case  the  difference  was  between  3 and  S%.  Results  are  based  on  two 
gages  at  the  610  and  700  metre  ground  ranges  and  upon  one  gage  at  the  794 
metre  range. 

PLUMBBOB,  Franklin-2.  For  this  shot  we  used  the  BRL  result  for  the 
single  ground  range  for  which  there  appeared  to  be  a good  measurement. 

PLUMBBOB,  Wilson-4.  For  this  shot  BRL  also  provided  values  of 
dynamic  pressure  impulse  as  listed  in  Table  2.  (We  integrated  the  waveform 
curves  as  a check,  however.) 

PLUMBBOB,  Hood-6.  Of  the  6 ranges  at  which  there  were  gage  results 
the  data  appeared  to  be  valid  at  4 locations.  There  were  two  gages  at  each  of 
these  locations.  We  used  the  BRL  results  as  listed  in  Table  2 (again  after  an 
integration  check). 

PLUMBBOB,  Kepler-9.  In  this  case  waveforms  obtained  at  three  ground 
ranges  were  quite  poor.  In  one  case  there  was  a large  disparity  between  two 
BRL  gage  results  (almost  a factor  of  4) ; in  a second  case  there  was  almost  a 
factor  of  2 difference  between  our  gage-reduced  result  and  BRL's  result.  At 
the  third  ground  range,  762  metres,  the  waveform  was  somewhat  better.  The 
percent  difference  between  our  result  and  BRL's  result  was  rather  large,  about 
20%;  in  Table  2 we  list  the  average  of  the  two  results.  We  discarded  the 
results  at  the  other  two  ground  ranges. 

PLUMBBOB,  Owens- 10.  For  this  shot  we  again  used  the  average  of  our 
gage-reduced  result  and  the  BRL  result  at  two  of  the  three  ground  ranges,  30S 
and  518  metres.  At  the  305  metre  range  the  results  are  in  reasonably  good 
agreement  (within  9%);  at  the  518  metre  range  the  two  results  differ  by  a 
factor  of  1.4  but  the  waveforms  are  quite  poor.  At  the  457  metre  range  our 
value  of  dynamic  pressure  impulse  was  exactly  equal  to  BRL’s  value  of  scaled 
dynamic  pressure  impulse;  the  BRL  value  appeared  anomalous  on  a data  plot, 
probably  as  a result  of  failure  to  scale  the  value  obtained.  In  this  case  we 
used  our  gage-reduced  result. 

2.2  Analysis  of  Tank  Data. 

We  next  consider  the  set  of  shots  listed  in  Section  1 on  which  tanks 
were  exposed  but  there  were  no  dynamic  pressure  measurements.  Our  rationale 
here  is  based  on  a finding  we  have  discussed  in  detail  in  Reference  1:  the 

displacement  which  a vehicle  exposed  to  a blast  wave  suffers  can  be  used  as 
a measure  of  the  dynamic  pressure  impulse  it  receives,  .'hat  is,  the  vehicle 
can  serve  as  a gage  for  dynamic  pressure  impulse  if  we  are  able  to  "calibrate" 
this  "gage".  The  calibration  is  the  curve  of  dynamic  pressure  impulse  versus 
displacement.  (We  are  ignoring  low  yield  devices  where  diffraction  effects 
play  a role.)  We  have  analyzed  the  procedure  and  results  for  1/4  ton  and  2 4 
ton  trucks  in  considerable  detail  in  Reference  1. 

We  would  not  be  able  to  use  the  procedure  for  tanks  if  we  had  no  dynamic 
pressure  impulse  data  for  the  ground  ranges  at  which  tanks  were  exposed.  That 
is,  we  would  have  no  way  of  obtaining  the  necessary  calibration  curve.  For- 
tunately, we  have  the  necessary  dynamic  pressure  impulse  data.  It  is  supplied 
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by  information  from  TEAPOT,  Met-12  and  Apple  11-13  and  from  UPSHOT-KNOTHOLE, 
Annie- 1 and  Grable-10.  The  data  for  the  pair  of  UPSHOT- KNOTHOLE  events  are 
gage  data  supplied  by  BRL;  for  Grable-10  we  also  have  results  based  on  both 
gage  and  (truck)  displacement  from  Reference  1,  but  at  larger  values  of  ground 
range  than  those  at  which  tanks  were  exposed.  The  BRL  data  and  the  Reference  1 
data  are  very  compatible,  i.e.,  very  well  represented  by  a single  curve  of 
scaled  dynamic  pressure  impulse  versus  scaled  ground  range  (as  shown  in  Figure 
7,  to  be  discussed  later).  The  results  are  listed  in  Table  2. 

For  Met-12  we  have  the  results  from  Reference  1,  listed  in  Table  2,  based 
on  both  gage  and  (truck)  displacement  data.  (There  is  a considerably  larger 
body  of  truck  displacement  data  than  of  tank  displacement  data  so  that  the 
dynamic  pressure  impulse  - displacement  curve,  i.e.,  the  calibration  curve,  is 
better  determined  for  trucks.) 

For  Apple  11-13  the  results  listed  in  Table  2 are  taken  only  partly  from 
Reference  1.  In  this  instance,  values  of  dynamic  pressure  impulse  at  the  518 
and  625  metre  ranges  listed  in  Reference  1 (based  only  upon  gage  data)  have 
been  averaged  with  BRL  gage  results,  i.e.,  the  two  results  given  equal  weight. 
At  the  625  metre  range  this  .makes  a negligible  difference  while  at  the  S18 
metre  range  the  value  of  dynamic  pressure  impulse  listed  in  Reference  1 differs 
from  that  of  Table  2 by  about  16%.  Finally  at  the  808  metre  range  the  result 
listed  in  Table  2 differs  from  that  listed  in  Reference  1;  the  latter  is  in 
error  as  a check  of  our  (previous)  data  reduction  reveals. 

2.2.1  Dynamic  Pressure  Impulse  Versus  Displacement  For  Tanks. 

We  plotted  dynamic  pressure  impulse  versus  displacement  for  the 
tanks  exposed  on  events  Met-12,  Apple  11-13,  Annie-1  and  Grable-10.  This 
provides  a total  of  14  data  points*  listed  in  Table  3 along  with  the  tank 
displacement  data  for  other  shots. 

There  are  several  variables  which  cause  deviations  of  the  data 
relative  to  a smooth  curve  through  the  points:  (I)  vehicle  orientation  with 

respect  to  the  bomb  — 4 data  points  correspond  to  side-on  orientation,  7 to 
face-on  orientation,  1 to  rear-on  orientation,  1 to  face-on  45®  orientation, 
and  1 to  face-on  3/4  left  orientation;  (2)  there  are  three  different  tanks  -- 
M4A3,  M24,  and  M48;  (3)  there  are  two  surface  conditions  --  rough  sand  and 
fine  sand.  There  obviously  are  not  sufficient  data  to  sort  out  the  effects 
of  the  variables.  However,  when  we  plotted  the  data  we  found  that:  (1) 

(initial)  orientation  appears  to  make  little  difference,  any  effect  being 
submerged  in  effects  of  the  other  variables;  (we  intuitively  expect  initial 
orientation  to  make  less  difference  for  tanks  than  for  trucks);  (2)  any 
systematic  deviation  which  could  be  attributed  to  type  of  vehicle  could  as 
well  be  attributed  to  inaccuracy  in  dynamic  pressure  impulse  --  for  example, 
the  Met-12  points  (points  B,  6,  and  17  in  Table  3)  are  a little  high  and  are 


v3 


* We  have  not  used  Point  8 of  Table  3,  a Smoky- 15  point.  Placing  this  point 
on  Figure  1 shows  that  there  is  clearly  something  wrong  with  it  --  and  what 
is  wrong  with  it  involves  the  displacement,  not  the  dynamic  pressure  impulse, 
even  though  the  latter  involves  extrapolation  on  Figure  6. 
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M48  points  while  the  Apple  11-13  points  are  a little  low  and  pertain  to  both 
M48's  and  K24's.  (In  fact,  if  we  use  an  average  dynamic  pressure  impulse 
versus  scaled  ground  range  curve  for  the  two  shots,  Met-12  and  Apple  11-13, 
these  deviations  disappear.) 

Plots  of  the  data  and  the  fits  obtained  are  shown  in  Figure  1.* 

We  used  these  data  to  infer  dynamic  pressure  impulse  from  the 
measured  displacements  on  events  GREENHOUSE  Easy-2;  UPSHOT-KNOTHOLE  Nancy-2, 
Badger-5,  and  Simon-7;  TUMBLER-SNAPPER  Fox-6  and  How-7. 

We  have  noted  in  Reference  1 that,  despite  rather  wide  fluctuations 
in  displacement  owing  to  the  several  variable  factors  mentioned,  the  central 
curve  following  the  trend  of  the  data  points  provides  a rather  good  calibration 
which  enables  us  to  use  a measured  displacement  to  determine  the  dynamic  pres- 
sure impulse  to  which  the  vehicle  was  subjected.  Since  displacements  vary  owing 
to  uncontrolled  factors,  the  inferred  dynamic  pressure  impulses  will  exhibit 
considerable  dispersion.  (However,  in  many  cases  described  in  Reference  1,  the 
displacement-inferred  values  are  about  as  reliable  as  the  gage- inf erred  values, 
especially  when  there  are  two  or  more  values  which  can  be  averaged.) 

The  calibration  curve  for  tanks  is  less  well  determined  than  that 
for  trucks  (Reference  1),  there  being  fewer  data  points.  Also  the  displacements 
are  smaller  resulting  in  greater  errors  especially  from  small  yield.  Nonetheless, 
when  the  displacement-inferred  dynamic  pressure  impulse  values  are  plotted  versus 
scaled  ground  range  (along  with  a few  data  points  for  which  there  also  are  gage 
data),  we  see  that  the  curve  obtained  is  reasonably  well  determined.  See  Figure 
3. 


We  used  an  eye-drawn  curve  rathar  than  a proportional  fit  in  Figure  1 
in  inferring  dynamic  pressure  impulse  from  measured  displacements  for  the  follow- 
ing reasons: 

(1)  point  20,  very  small  displacement.  Is  much  less  important 
than  the  other  data  points  in  application  of  this  fit,  i.e.,  inference  of 
dynamic  pressure  impulse  from  measured  displacement; 

(2)  point  21  has  virtually  the  same  displacement  as  the  cen- 
troid so  that  it  can  be  ignored  in  drawing  a straight  line  through  the  data 
points;  [the  line  must  pass  through  the  centroid: 

Zfnr  - 2.448,  In  1 = 3.768; 

q 

on  the  abscissa  and  ordinate  scales  from  Figure  1 we  see  that  both  fits  do 
pass  through  the  poiht  5^5-55  . u 56/  IqIErT^-  . 43.29;  the  line,  do  not 

pass  through  the  poiht  ? . 21-85  - . 

q 

* Because  of  the  mentioned  variables  whose  effects  we  cannot  disentangle  owing 
to  the  small  amount  of  data,  we  have  numbered  the  points  on  Figure  1.  Cor- 
responding numbers  are  listed  in  Table  3.  By  comparison  one  can  verify 
that  the  data  do  not  exhibit  any  marked  effects  owing  to  differences  among 
the  above  variables. 
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(3)  the  Grable-10  points  21  and  22  are  badly  inconsistent  with 
one  another  and  points  9 and  10  are  somewhat  inconsistent.  Thus  our  eye-drawn 
curve  was  a little  higher  than  Fit  1 of  Figure  1 at  the  small  displacement  end 
and,  like  Fit  1,  passed  through  the  centroid  of  the  data  points.  The  differ- 
ence, however,  is  nowhere  greater  than  a few  percent.  We  also  note  from 
Figure  1 that  a curve  through  the  data  points  would  do  no  better  than  a 
straight  line  and  is  therefore  not  warranted.  It  is  quite  possible,  however, 
that  if  the  data  covered  a wider  range  a straight  line  would  no  longer  be 
adequate. 

2.3  Discussion  Of  and  Minor  Improvements  Upon  Previously  Used  Blast  Data. 

Most  of  the  dynamic  pressure  impulse  data  listed  in  Table  2 which  we 
have  not  yet  discussed  is  taken  from  Reference  1.  In  a few  instances  a value 
taken  from  Reference  1 (gage  only  - no  vehicle  displacement  data  available) 
was  averaged  with  a BRL  gage-reduced  result  as  shown  in  the  table.  In  one 
instance,  Turk-4,  ground  range  595  metres,  we  reduced  the  gage  data  even 
though  the  gage  had  cut  off  near  the  end.  We  estimated  the  shape  and  time  of 
pulse  completion  from  other  Turk-4  data  and  believe  that  the  inaccuracy  so 
introduced  is  small  — not  more  than  a few  percent.  Reference  to  Figure  3 
shows  why  this  is  desirable:  Turk-4  has  the  lowest  scaled  burst  height  for 

which  we  have  gage  data  (except  for  a single  Annie- 1 point  and  except  for 
surface  burst  data);  its  scaled  burst  height  is  only  moderately  greater  than 
those  of  other  events  shown  on  Figure  3;  the  Turk-4  data  are  consistent  with 
the  data  for  the  other  events  shown  on  Figure  3 and  extend  to  a much  greater 
value  of  scaled  ground  range.  Thus  the  curve  fit  of  Figure  3 is  not  diminished 
in  reliability  by  inclusion  of  Turk-4  data  which,  however,  enable  its  use  over 
a much  greater  range  in  constructing  height  of  burst  charts,  the  object  of  this 
report. 

For  events  Yuma-4  and  Wasp  Prime-9  we  used  our  values  from  Reference  1, 
therein  described  as  "first  iteration"  results;  i.e. , when  available,  truck 
displacement- inferred  dynamic  pressure  impulse  values  were  averaged  with  gage 
values  in  obtaining  our  best  estimates  of  dynamic  pressure  impulse.  (The 
manner  of  averaging,  justification,  and  discussion  of  accuracy  and  reliability 
are  considered  in  detail  in  Reference  1.)  We  originally  considered  Yuma-4  to 
be  an  ideal  event,  but  we  now  believe  the  values  achieved  by  averaging  displace- 
ment-inferred dynamic  pressure  impulse  values  with  the  gage  values  represent  an 
improvement  over  the  gage  values  alone;  the  results  achieved  are  certainly  more 
compatible  with  the  other  data  in  the  same  scaled  height  of  burst  region.  See 
Figure  9 (and  compare  with  gage  only  data  listed  in  Reference  1,  page  115). 

For  event  Encore-9,  a near- ideal  event,  we  use  the  gage  results  --  even 
though  in  Reference  1 we  sought  to  improve  upon  these  results  by  use  of 
displacement -inferred  dynamic  pressure  impulse.  The  reasoning  used  in  Reference 
1 was  that  since  the  data  scatter  on  a displacement  versus  dynamic  pressure 
impulse  plot  was  just  as  great  for  idea  1/near- ideal  shots  as  for  non-ideal  and 
since  the  displacement- inferred  dynamic  pressure  impulse  for  non- ideal  shots 
is,  on  the  average,  as  accurate  and  reliable  as  the  gage  values  --  then 
averaging  of  displacement- inferred  with  gage- inferred  dynamic  pressure  impulse 
should  improve  the  values  for  near- ideal  events  just  as  it  does  for  non- ideal. 
The  points  at  issue  are  discussed  in  detail  in  Reference  1.  One  point,  however, 
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is  that  there  is  a much  greater  variety  of  ideal/near- ideal  shots  than  of 
non- ideal,  especially  regarding  the  range  of  weapon  yields  in  the  data  base. 
When  good  gage  data  are  available,  these  data  should  produce  better  results 
than  can  be  obtained  by  averaging  in  displacement- inferred  values  --  owing  to 
the  effects  of  the  several  uncontrolled  variables  previously  described.  That 
is,  while,  on  the  average,  displacement-inferred  values  are  reliable,  consider 
able  dispersion  is  to  be  expected.  Compare,  for  example.  Figure  3;  here  with 
the  exception  of  the  Turk-4  data  points  and  a single  Annie- 1 point,  all  of  the 
plotted  points  are  based  entirely  upon  displacement- inferred  values  of  dynamic 
pressure  impulse.  While  the  central  curve  following  the  general  trend  of  the 
data  is  reasonably  well  determined,  the  data  point  deviations  are  rather  large 
— considerably  larger  than  the  deviations  in  Figures  4-10. 

Near- ideal  shots  also  have  (small)  real  differences  relative  to  one 
another.  This  is  shown,  for  example,  by  Figure  3.13  of  Reference  1 which  is 
a plot  of  scaled  dynamic  pressure  impulse  versus  scaled  ground  range.  Data 
for  the  near- ideal  (surface  burst)  shots  do  not  completely  coalesce  under  the 
scaling  --  there  are  clearly  small  but  real  differences  owing  to  variable 
factors  not  controlled  in  the  experiments.  In  Figures  3-10  the  data  for  the 
events  plotted  on  each  figure  seem  to  coalesce  quite  well  with  respect  to  a 
single  curve.  In  Figure  7 some  deviation  can  be  seen.  In  general,  deviations 
tend  to  appear  when  data  for  several  events  are  shown  on  a single  plot,  when 
there  are  several  data  points  for  each  event,  and  where  the  data  for  several 
events  overlap,  i.e.,  cover  the  same  domain  of  the  abscissa,  scaled  ground 
range.  The  data  in  Figures  3-10  pertain  to  a much  smaller  range  of  weapon 
yields  than  do  the  surface  burst  data  of  Reference  1 (Figures  3.13  and  3. IS). 
In  general,  when  the  data  for  sever;  1 shots  coalesce  we  can  attach  a high 
degree  of  reliability  to  the  data. 

In  constructing  height  of  burst  charts  in  this  report  we  used  the  results 
described  and  listed  in  Table  3.  For  the  surface  burst  data  we  used  Figure 
3.13  of  Reference  1.* 

Finally  we  analyzed  field  data  for  many  more  shots,  but  found  the  data 
unusable.  In  some  instances  the  waveforms  were  very  bad;  in  other  instances 
the  gages  were  placed  at  elevations  other  than  3 feet  above  the  ground.  The 
effect  of  gage  elevation  is  non-negligible.  (Different  types  of  gages  also 
exhibit  somewhat  different  responses  which  must  be  accounted  for  in  order  to 
achieve  consistent  results.)  - 

NOTE:  The  BRL  data  used  in  Table  3 have  not  been  published. 

2.4  Dynamic  Pressure  Impulse  Versus  Displacement  for  Self-Propelled  Howitzers 

Figure  2 is  a plot  of  displacement  - dynamic  pressure  impulse  data  for 

* For  a single  (extrapolated)  point  on  the  charts  for  which  the  scaled  dynamic 
pressure  impulse  is  15  kPa-sec  and  the  scaled  height  of  burst  ’is  zero,  we 
also  used  Figure  3.15  of  Reference  1 and  took  account  of  the  difference 
between  Figures  3.13  and  3.15  in  this  region.  Figure  3.15  extends  to  slight 
ly  higher  values  of  dynamic  pressure  impulse  than  does  Figure  3.13  so  that 
less  extrapolation  is  required. 
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self-propelled  howitzers.  Here  we  used  values  of  dynamic  pressure  inpulse 
already  described.  There  are  only  eight  data  points.  These  involve  2 vehicles. 
3 orientations,  2 soil  types  and  7 nuclear  events.  Three  of  the  data  points 
(2  Badger-5  and  1 Siaon-7)  result  from  shots  for  which  there  were  no  dynamic 
pressure  neasurenents.  Therefore,  the  values  we  used,  taken  from  Figure  3 and 
from  the  scaled  height  of  burst  chart,  Figure  11,  are  actually  based  upon  the 
displacement- inferred  dynamic  pressure  impulse  values  for  tanks. 

We  did  not  attempt  to  use  the  howitzer  displacement  data  to  infer  values 
of  dynamic  pressure  impulse.  To  do  so,  we  should  proceed  in  the  sane  manner 
as  for  the  tank  data.  Referring  to  Figure  2,  we  see  that  of  the  eight  data 
points  only  five  (points  31-35)  could  be  used.  No  curve  of  any  reasonable 
degree  of  reliability  can  be  drawn  based  on  these  points. 
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SECTION  3 

(This  Section  is  Unclassified) 

RESULTS  AND  CONCLUSIONS 

The  results  described  are  given  in  Tables  2 and  3.  Least  squares  fits 
are  given  in  Tables  4 and  5.  A description  of  the  methods  used  in  obtaining 
these  fits  and  the  definitions  and  explanations  of  the  several  measures  of 
goodness  of  fit  listed  are  given  in  the  Appendix.  The  simplest  measure  of 
goodness  of  fit  to  understand  is  the  root  mean  square  percent  error  (E^  and 
Ej  ).  From  Table  4 we  see  that  for  the  displacement  - dynamic  pressure  impulse 
q fits,  the  errors  are  quite  large.*  The  error  in  estimating  dynamic  pres- 
sure impulse  from  given  displacement  is  fairly  large  (35-40%)  but  tolerable 
and  as  we  have  seen  the  values  obtained,  on  the  average,  are  reasonable  and  of 
considerable  value  especially  when  gage  data  are  poor  or  lacking.  If  we  were, 
however,  to  attempt  to  infer  displacement  from  a given  dynamic  pressure  impulse 
the  errors  would  be  much  larger  (^150%)  — though  again  the  results  would  repre- 
sent reasonable  approximations  to  the  average  displacement  to  be  expected.  (The 
reason  for  the  wide  disparity  between  the  errors  for  the  two  types  of  inference 
is  seen  to  be  simply  due  to  the  slope  of  the  fitted  lines  on  the  log-log  plots. 
Figures  1 and  2,  or  equivalently  from  the  exponents  in  the  fits  listed  on  the 
Figures  and  in  Table  4.) 

The  percent  errors,  Ej  , shown,  in  Table  5 for  the  dynamic  pressure  impulse 
versus  scaled  ground  range  ^ fits  are  quite  modest,  varying  from  about  4 to  25%. 
These  errors  are  quite  comparable  to  those  obtained  in  Reference  1 for  the  var- 
ious fits  therein. 

The  results  described  are  shown  along  with  the  fits,  in  Figures  1 - 10. 
These  results  were  then  used,  along  with  the  surface  burst  data  in  Reference  1 
(Figure  3.13  and  in  one  instance  Figure  3.15)  to  construct  the  desired  scaled 
dynamic  pressure  impulse  contours  as  shown  in  Figures  11  and  12.  Also  shown 
is  the  locus  of  points  separating  the  regular  and  Mach  reflection  regions. 

We  have  drawn  the  contours  to  conform  as  accurately  as  possible  to  the 
plotted  points  while  maintaining  smoothness  and  a continuous  variation  in 
contour  shape  as  we  proceed  from  low  to  high  values  of  scaled  dynamic  pressure 
impulse. 

Figure  11  is  the  better  and  more  reliable  of  the  two  charts.  The  curves 
fit  the  data  points  much  better  than  for  Figure  12.  In  the  latter  case  we  have 
not  attempted  to  force  the  curves  through  the  data  points  as  this  would  lead  to 
structure  in  the  contour  shapes  which  is  not  justified  by  the  data. 

There  are  various  degrees  of  dust**,  the  Met-12  data  corresponding  to  the 
heaviest  dust  case  while  the  Apple  11-13  and  Bee-6  data  represent  more  moderate 


* The  errors  here  are  much  larger  than  in  Reference  1 because  there  are  much 
more  data  for  trucks  (of  concern  in  Reference  1)  than  for  the  tanks  and 
howitzers  of  concern  here. 

**  We  now  believe  that  the  degree  of  "dustiness"  is  an  indicator  of  the  severity 
of  the  precursor  which  is  related  in  turn  to  the  temperature  of  the  pre- 
shock thermal  layer, 
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dust.  This  is  the  principal  cause  of  the  difficulty  in  Figure  12.  However,  we 
do  not  believe  that  separate  charts  for  moderate  dust  and  for  heavy  dust  are 
warranted.  This  distinction  is  simply  too  fine  considering  the  quality  of  the 
data. 


In  conforming  the  contours  to  the  data  points  in  as  reasonable  a manner 
as  possible,  we  note  that  the  two  charts  differ  for  the  higher  scaled  burst 
heights,  i e.,  the  comparable  contours  do  not  coincide,  even  though  the  same 
data  are  used  for  the  two  highest  scaled  burst  heights  on  Figures  11  and  12. 

We  cannot  regard  the  contours  in  this  region  as  well  determined;  our  preference 
in  this  region  is  for  Figure  11,  because  of  the  greater  smoothness  of  the  con- 
tours and  because  we  believe  the  data  at  the  third  highest  scaled  burst  height 
are  better  for  Figure  11  than  for  Figure  12. 

Note:  In  applications  of  Figures  11  and  12  the  scaling  factor  S,  (see 

General  Nomenclature)  should  be  used  rather  than  W“^/3  for  scaling  ground  range 
and  height  of  burst;  Si  should  be  used  in  conjunction  with  contours  rather  than 
IK-1/3  for  instances  in  which  ambient  pressure  and  temperature  are  specified  and 
differ  significantly  from  the  standard  values  of  14.7  psi  (101.4  kPa)  and  15°C, 
respectively. 
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SECTION  4 

(This  Section  is  Unclassified) 

REOWtENDATION 

It  is  reconaanded  that  the  Height  of  Burst  Chart,  Figure  11,  be  accepted 
as  standard,  i.e. , as  the  acst  accurate  obtainable  from  the  totality  of 
existent  blast  wave  measurements.  Figure  12  which  includes  Iq  for  moderate/ 
heavy  dust  environments  is  for  unique  conditions.  The  description  of  these 
conditions  was  discussed  in  Reference  1.  Thus,  Figure  12  should  be  used  only 
for  siailar  conditions  or  be  ignored. 
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APPENDIX  A 

(This  Appendix  is  Unclassified) 
THE  LEAST  SQUARES  PROCEDURES 

1.  INTRODUCTION. 


All  of  the  least  squares  fits  in  this  report  were  obtained  with  the  TI-59 
calculator  using  library  programs  supplied  by  Texas  Instruments. 

All  of  the  fits  were  either  of  the  form 

y ■ ax  ♦ b (1) 

or  y • aQ  ♦ a1x1  ♦ a^  (2) 

where  y was  the  natural  logarithm  of  the  desired  quantity;  x was  either  the 
independent  (controlled)  variable  or  its  logarithm  and  the  least  squares 
variables  Xj,  x2  involved  only  the  special  case  x1  * x and  x2  * x2. 

The  least  squares  equations  are  obtained  by  minimising  the  sum  of  the 
squares  of  the  deviations  of  the  data  from  the  curve.  If  the  experimental 
points  have  a variable  scatter  for  a given  small  range  of  x,  weights  should 
be  applied  in  the  procedures  so  that  the  least  squares  equations  correspond- 
ing  to  Equation  (1),  for  example,  are  obtained  by  setting  the  derivatives  of 

S • ZWjC/j  - (mxj  ♦ b)32  C3) 


with  respect  to  m and  b equal  to  zero,  thus  obtaining  two  linear  algebraic 
equations  which  can  be  solved  for  a and  b.  Here  S is  the  weighted  sum  of  the 
squares  of  the  residuals  about  the  fitted  curve  and 


Wi  * °2/oi2 

where  a is  a constant  (to  be  determined  from  the  deviations  of  the  data  from 
the  fitted  curve)  and  o^  is  a measure  of  the  expected  deviation  from  the 
true  value  for  an  observation  y^  (of  unit  weight). 

In  the  Texas  Instruments  programs  ■ 1.  In  individual  cases  to  be 
discussed  we  will  point  out  that  the  progress  used  are  nonetheless  quite 
adequate  for  our  needs.  This  is  mainly  due  to  the  fact  that,  in  all  cases, 
we  minimized  the  sum  of  the  squares  of  the  logarithm  of  the  desired  function. 
This  is  equivalent  to  minimizing  the  sum  ot  the  squares  of  the  percentage 
deviations  from  the  curve  since 


A(6t  y)  ■ ^ 


(4) 
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That  is,  for  small  deviations  from  the  cut/e,  the  fractional  error  & is 
equal  to  A(6i  y) . For  large  deviations  this  is  n.»t  exact  and  in  our  tabula 

tion  of  results  we  show  both  , 

*'  y “ y ** 

£(6t  yi  • In  y -ri  — ) 

i ; } c 


where  yc  ■*  the  value  of  yfx^)  corresponding  to  the  curve  fit.* 

The  statistical  analog  of  Equation  (4)  (Reference  3)  is 

Var  (In  y)  ■ -7.  Var  y (Var  ■ variance)  . (4a) 

y 

Now  in  all  of  the  cases  the  values  of  y covered  a large  range.  If  we  con- 
sider, for  example,  a displacement  measurement  of  100  metres  to  be  in  error 
by  2%,  the  square  of  its  deviation  from  the  true  value  is  4;  if  we  consider  a 
displacement  measurement  of  1 metre  to  be  in  error  by  20%,  the  square  of  its 
deviation  from  the  true  value  is  0.04.  So  even  though  the  percent  errors  in 
our  data  are  not  necessarily  uniform  (as  best  these  errors  are  known) , assum- 
ing thorn  to  be  uniform  is  much  closer  to  reality  than  any  other  assumption  we 
can  make.  From  the  example  just  given  we  readily  see  that  the  coefficients 
in  the  least  squares  fits  would  be  determined  almost  entirely  by  the  data 
with  large  values  of  y if  we  were  to  minimize  absolute  rather  than  percent 
deviations  from  the  fitted  curves. 


2.  DYNAMIC  PRESSURE  IMPULSE  VERSUS  DISPLACEMENT  (a) 

DISPLACEMENT  VERSUS  DYNAMIC  PRESSURE  IMPULSE  (b) 

In  application  of  our  results  we  need  fits  of  the  data  for  cases  (a)  and 
(b),  i.e.,  with  each  quantity  used  in  the  role  of  independent  (controlled) 
and  of  dependent  variable. 

In  each  subcase  of  this  case  the  data  are  fitted  very  well  by  straight 
lines  on  log-log  plots.  Thus  for  fitting  on  the  TI-59  calculator  the  fit  is 
of  the  form  y ■ mx  ♦ b (1) 

with  y ■ in  Iq,  x ■ in  D for  Case  (a) 

and  y ■ in  D,  x ■ in  I-  for  Case  (b) 


and 


Lq' 

in  D, 
b » in  B 


(S) 


* In  most  cases  we  omit  the  subscript  i denoting  the  ith  data  point;  the 
summations  are  taken  over  the  data  points  in  all  cases  so  that  omission 
of  the  subscript  1 will  not  lead  to  any  confusion. 
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Equation  (1)  can  then  also  be  written  [for  Case  (a)],  1^  • BDm  (la) 

with  a similar  relation  for  Case  (lb).  However,  after  we  have  fitted  Case  (a) 

we  can  calculate  (b)  directly  (without  fitting  it)  provided  we  recall  certain 

data  summations  from  the  machine  memory.  For  Equation  (1)  we  obtain  from  the 

TI-59  fit:  r r r /u 

slope  - m - $Xl 


Exfc  - (Ex)‘/N 


Yl  . b - fr  ~ aE* 

7 intercept  n 

, 1/2 


(6) 


(7) 


and  the  correlation  coefficient,  r ■ m (-£*■)  (8) 

ay 

r2  is  called  the  coefficient  of  determination.  The  various  summations  can  be 

2 2 

recalled  from  the  machine  memory.  ax  and  oy  are  the  variances  of  the 
x- array  and  y- array  data  and  are  given  by  * 


°x  ■ if  ^ 

(9) 

■ ff  - ?>2 

(10) 

with 

x • ^ Ex 

(11) 

and 

7 ■ s cy 

(12) 

N is  the  number  of  data  points. 

Equations  (6)  to  (8)  are  not  symmetrical  in  x and  y and  the  algebraic 
inverse  of  the  least  squares  fit  I 


n 


is 


0 . <i) 


BD 


1/m  1/m 

I 


(la) 


but  this  is  not  a least 


squares  fit  to  the  data  with  the  roles  of  Iq  and  D reversed. 

However,  manipulation  of  Equations  (6)  to  (10)  leads  to  the  equation 

r - [Exy  - I ExEy]  ([Ex2  - ± (Exl2][Ey2  - i (Ey)2]]  U 2 (3a) 

which  is  symmetrical  in  x and  y,  i.e.,  unchan'-  when  x and  y are  interchanged. 


* Bar  over  a quantity  indicates  an  average  value. 
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Thus,  if  we  have  obtained  the  least  squares  fit  Iq 


* BDm  we  know  that 


the  correlation  ratio  r has  the  same  value  with  the  roles  of  Iq  and  D 
reversed.  It  can  be.  shown  (by  taking  the  origin  at  the  data  centroid) , 
(References  3 and  4),  that  the  least  squares  fits 


BD 


AI 


.m 


t m ■ r‘ 


(la) 

(lb) 
(13) 


have  values  of  l and  in  which  satisfy 

Thus,  having  determined  m (and  r)  in  the  least  squares  fit  to  case  la 
we  obtain  l in  the  least  squares  fit  to  Case  (lb)  from  Equation  (13).  The 
constant  A in  Equation  lb  is  then  determined  by  the  fact  that  both  least 
squares  fits  (la)  and  (lb)  pass  through  the  centroid  of  the  data  distribu- 
tion. (Compare  Reference  S.)  Thus,  from  Equation  (lb)  we  have 

in  V • in  L in 

and  in  particular  tn  D • in  A+£  in  1^ 

where  from  Equations  (11)  and  (12) 

Tn~0  - £ Utn  D) 

Iq) 

and  these  values  are  available  from  the  machine  memory  after  running  case 
(la).  Thus  from  Equation  (14)  we  determine  In  A,  hence  A and  the  least 

squares  fit  to  case  (lb)  is  fully  determined. 

Having  obtained  the  least  squares  fits,  the  coefficient  of  determination 
provides  a measure  of  the  goodness  of  fit.  By  manipulation  of  Equation  (3), 
with  ■ 1,  and  Equations  (6)  through  (12)  it  can  be  shown  that 


(lb) 

(14) 

(15) 


(16) 


NO. 


In  this  report  we  are  concerned  only  with  positive  correlation,  so  that  the 
nearer  r is  to  *1.  the  smaller  is  S,  the  sum  of  the  squares  of  the  deviations 
from  the  fitted  curve,  and  hence  the  better  the  fit  to  the  data.  The  expres- 
sion for  1 - r2  can  also  be  put  in  the  form 
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l - r - — — 

Uy  - y)‘ 

thus  expressing  i - r2  as  a ratio  of  the  sum  of  the  squares  of  the  deviations 
from  the  curve  to  the  sum  of  the  squares  of  the  deviations  from  the  mean  (of 
the  data).  The  denominator  may  be  regarded  as  a normalizing  factor.  It  s 
independent  of  the  functional  form  used  for  the  fitting  function.  Without 
such  a denominator  1 - r2  would  tend  to  increase  with  the  number  of  data  points 
— even  if  the  data  were  excellent. 

Now  letting  J and  D be  the  sum  of  the  squares  of  the  deviations 

from  the  fitted  curve  when  in  Iq  and  in  D,  in  turn,  play  the  role  of  dependent 
variable,  we  have  from  Equation  (16) 

(i  - r2)  - S£n  *q  - ° . (18) 

NoL  Iq  No£n  D 

Here  ain  I » ffttC*  V2  * ^ V2] 

q (19) 

<4  d - £ esc*  °)2  - jf  CS*  D)2] 

in  accord  with  Equations  (9)  through  (12).  Since  we  have  seen  that  r is 
unchanged  when  the  roles  of  D and  Iq  ere  reversed.  Equation  (18)  shows  that 
the  sums  of  the  squares  of  the  deviations  from  the  curves  are  not  the  same 
for  the  two  corresponding  least  squares  fits  (to  the  same  data)  but  that 


Sin  D 


0 in  I 

q 

ain  D 


In  all  cases  treated  in  this  report  we  found  S^t  jq 


(20) 


to  be  considerably 


less  than 

In  an  application  to  be  made  of  the  results  of  this  case  we  wish  to 
know  the  standard  deviation  of  the  value  of  the  slope  of  the  curve  as  well  as 
the  standard  deviation  of  our  observation.  We  use  the  following  notation 
with  the  subscript  c referring  to  the  curve  fit  in  each  instance: 
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at  T - • variance  of  an  observation  of  in  Ia  relative  to  the 
cn  ig i c * 

M curve;  (note  that  we  have  assumed  the  variance  of 

the  percent  error  to  be  independent  of  the  abscissa 

value) 

-2  2 _ . _ r-  _ 1 _ ... 


,,  0“  * variance  of  slopes  of  curve  fits, 

c m 


Then  it  can  be  shown  that  (Reference  6) 


^ Iq.C  * s 


/ 

^ 'J 


hn  D ,C  * *ln  0/ 

2 W"2 


0/  * — 

* NO,. 


1 ‘ r °fn  D 


/ N‘ 
ia/ 


1 ‘ 1,2  °f  n I 


r" 


Using  Equation(8)  twice,  first  with  D,  I being  independent  and  dependent 
variable  and  then  with  their  roles  reversed,  we  also  see  that 

a m ' 

m 

In  the  tabular  results  for  the  Displacement  - Dynamic  Pressure  Impulse 
data  curve  fits,  we  list  for  each  subcase  the  quantities  N,  j * i„,c 

S.  rv.  a,  t,  m,  o.,  a in  addition  to  r and  the  fitting  functions. 

Ot  u cti  d,u  c m 

To  gain  further  insight  into  the  reliability  of  the  data  fits  we  have 
also  tabulated  the  following  additional  quantities  defined  as  follows: 


Here  there  are  N-2  degrees  of  freedom;  the  two  degrees  of  freedom  lost 
correspond  to  the  number  of  regression  coefficients  (two) . 
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0 - D, 


I.  - I, 


V * <-*£-> 


and  finally  the  root  mean  square  percent  error  which  is 


(S£/N)1/2  x 100% 


E7  - (Si  m x 100% 

q xq- 


The  o'  quantities,  Equation  (25),  are  standard  deviations  of  the  fractional 
error  (percent  error  apart  froa  a factor. of  100).  Thus  when  the  number  of 
data  points  N is  fairly  large  (large  enough  so  that  the  fit  is  fairly  reli- 


able) the  values  of  o'  oi  c are  almost  the  same  as  the  values  of  E-.  and 


^Iq'  respectively,  (aside  from  the  factor  of  100  expressing  the  latter  quan- 


tities in  percent) . 

The  reason  we  have  chosen  to  list  these  additional  quantities,  in  this 
and  in  other  data  fits  to  be  discussed  shortly,  is  that  in  many  cases  the 
data  scatter  is  rather  large  (irrespective  of  the  functional  form  selected 
for  fitting  the  data) . This  means  that  a deviation  y in  In  y for  a given 
data  point  may  differ  considerably  from  Ay/y.  To  illustrate,  suppose  the 
curve  fit  and  data  point  values  of  a displacement  are  10  metres  and  6 metres. 


The  contribution  of  this  point  to  the  sum  of  the  squares  of  the  deviations 

_ r f ..  k « . t .n  o _ a 1 - _ c . : . ..  1 


of  In  D is  then  (In  10-  In  6)  « 0.261  (irrespective  of  which  value  is  the 
curve  fit  value  and  which  is  the  data  point  value) . 


The  contribution  to  S 1 


(-^)  * 0-160 


if  the  curve  fit  value  is  10  and  the  data  point  value  is  6;  if,  however,  the 
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curve  fit  valut  is  6 and  the  data  point  value  is  10,  the  contribution  to 

0.444. 


is 


Since  the  least  squares  curve  tends  to  pass  through  the  data  region  with 
a fairly  uniform  distribution  of  points  on  either  side,  on  the  average  the 
curve  fit  will  be  below  the  data  point  about  as  often  as  above  it.  Averaging 

0.302  which  is  to  be  compared  with 


...  0. 160  ♦ 0.444 

the  above  results  yields  5 


the  value  0.261  for  the  logarithmic  deviation.  Had  we  chosen  the  values  10 
and  9.8  instead  of  10  and  6,  the  two  compared  numbers  would  be  virtually 
identical. 

It  is  inherent  in  the  nature  of  the  data  we  are  dealing  with  in  this 
report  that  there  is  large  data  scatter  while  the  number  of  data  points  is 
quite  small.  We  found  (Reference  1),  using  the  same  least  square  procedures 
described  herein,  that,  when  the  scatter  of  the  data  is  small)  i fairly 

close  in  value  to  Q and  similarly  SJ  is  fair-y  close  to  ^ . Usually 


S„  _ < SI  and  5.  T < SI  . The  situation  is  sometimes  reversed,*Siowever, 
•Ol  u D ZJt  I q i q 

usually  because  one  or  two  data  points  with  large  deviations  from  the  curve 


fit  are  below  the  curve.  Conversely,  for  a case  in  which  one  or  two  data 
points  are  far  above  the  curv*  it  can  happen  that  S^,  » y since  there  is 


no  limit  to  the  contribution  such  a point  can  make  to  S^,  y here  indicating 


any  quantity  whose  In  has  been  fitted. 


the  maximum  contribution  to  S’  is 


y-ye 


(For  a data  point  below  the  curve  fit 
2 


( -) 


o-y 


1 0 


The  results  shown  in  Table  4 are  consistent  with  these  statements  taken  from 
Reference  1.  [The  large  effect  which  a single  data  point  can  have  may  be 
seen  by  calculating  D-D 


(-g-£)  for  point  31  of  Figure  2.  The  value,  18.52,  is 
c 


almost  as  large  as  S£  (»  19.06)  for  the  8 data  points  of  Figure  2.] 


Finally,  quantities  such  as  oj  c and  Eiq  are  directly  associated  with 
the  plotted  data  and  are  easily  visualized.  Since  they  pertain  to  fractional 
or  percent  errors  they  apply  equally  well  at  all  points  along  the  curve  fit 
(although  the  curve  fit  has  greater  predictive  accuracy  in  the  vicinity  of 
the  data  centroid  than  toward  the  ends  of  the  data  range) . On  the  other 
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hand  £ refers  to  the  natural  logarithm  of  the  plotted  data  (ordinate) 

and  the  latter  may  be  quite  small  or  large  (or  even  negative),  so  that  while 
atn  Iq  a standard  measure  of  the  data  scatter  with  respect  to  the  fitted 
curve,  it  is  not  as  simple  to  interpret  as  are  the  additional  quantities 
which  we  have  tabulated. 

In  sum,  the  tabulated  measures  of  goodness  of  fit  provide  some  insight 
into  the  assessment  of  the  reliability  and  utility  of  the  data  and  the 
curve  fits  and  help  with  the  interpretation  of  results. 

The  quantities  described  here  are  also  tabulated  for  various  curve  fits 
of  data  other  than  displacement-dynamic  pressure  impulse  data.  We  discuss 
the  least  squares  treatment  of  the  remaining  data  in  the  following  section. 


3.  ADDITIONAL  LEAST  SQUARES  FITS. 

Scaled  Dynamic  Pressure  Impulse  versus  Scaled  Ground  Range: 

The  fitted  curve  is  given  by 

y • »q  ♦ ♦ a2x2  (27) 

where  x ■ scaled  ground  range,  y ■ In  S^I  . 

In  running  this  case  on  the  TI-S9  calculator  we  feed  in  the  natural 
logarithm  of  each  data  point.  Again,  the  machine  generates  a fit  and  all  of 
the  relevant  summations  over  the  data  are  available  from  the  machine  memory. 
For  this  type  of  fit  (trivariate) , however,  the  machine  provides  a quantity 
R2  rather  than  r as  a measure  of  the  goodness  of  fit.  R is  the  multiple 
linear  correlation  coefficient  between  y and  the  other  least  squares  variables 

2 2 

x,  ■ x and  x2  ■ x ; R is  called  the  coefficient  of  determination  (see  Ref- 
erence 6)  and  is  given  by 


Tyl  * ry2  ~ 2r12  ryl  ry2 


1 - r 


12 


(28) 


The  subscripts  1 and  2 refer  to  Xj  and  x2;  r j,  ry2,  and  r12  are  given  by 
equations  analogous  to  Equation  (8a) . 


yi 


[Exy  - £££]  ([Ex2  - -^^][Ey2 


(29) 
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ry2  “ 


[Ex2y 


#3  (ctx4  - ^][£y2  - V') 

12 . tfc*  . Sfi]  (cr,2  - iifliKix4  - fl^]) 


-1/2 


-1/2 


(29) 


(The  superficial  lack  of  symmetry  between  and  r^?  results  from  the  fact 
that  the  least  squares  variables  Xj  and  x2  as  used  here  are  x and  x2,  respec- 
tively.) The  quantities  r^,  ry2,  and  Tp  are  the  coefficients  of  correlation 
(also  called  staple  correlation  or  tero-order  coefficients)  between  y and  Xp 
y and  x2,  and  and  x2,  respectively.  The  important  point  here  is  that 


2 ztyycr 

(l  - r2)  V 

Z(y-y) 


(17a) 


so  that  R*  provides  a aeasure  of  the  goodness  of  fit  precisely  similar  to  that 
provided  by  r2  in  the  previous  discussion:  R2  is  a measure  of  the  closeness 

of  fit  of  the  regression  plane  (in  In  y,  Xj,  x2  space)*  to  the  data  points. 

There  are  two  further  changes:  replacement  of  N-2  by  N-3  in  Equations 

(21)  and  (25)  and  replacement  of  Iq  by 


S.I  ,C  * Sfn  S.I  /N’ 
i q i q/ 


aU  D,C  " SZn  u/H~: 


(21a) 


and 


°D,C  " ^S0/N-3> 


1/2 


JS  T C m CSS  I /N‘3> 

Vq’L  Vq 


1/2 


(2Sa) 


q'~  i q 

Referring  to  the  results  listed  in  Table  5 we  see  that  usually 
SI  < S’  («d  in  *lct  this  relation  holds  for  *11  of  the  cases  where 

i q siiq 

the  number  of  data  points  ^11).  The  two  values  are  usually  quite  close, 
however,  as  might  be  expected,  since  the  data  dispersion  is  fairly  modest. 


* *1  * x ■ scaled  ground  range  x2  ■ x‘ 
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TABLE  3.  (U)  SUMMARY  OF  DISPLACEMENT  DATA  FOR  TANKS  AND  HOWITZERS 
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TABLE  5.  (U)  LEAST  SQUARES  FIT  RESULTS 

Scaled  Dynamic  Pressure  Impulse  vs.  Scaled  Ground  Range 
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